Interaction-induced beats of Friedel oscillations in quantum wires 
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We analyze the spectrum of electron density oscillations in an interacting one- dimensional electron 
system with an impurity. The system's inhomogeneity is characterized by different values of Fermi 
wave vectors kp — k L i R on left/right side of the scatterer, leading to a Landauer dipole formation. 
We demonstrate, that while in the noninteracting system the Friedel oscillations possess only one 
periodicity related to the local kp, say kh on the left side, the interplay of the interactions and 
the Landauer dipole generates an additional peak in the spectrum of density oscillations at the 
counterpart kn. Being only present in correlated systems, the position and shape of this spectral 
feature, which in coordinate space is observable as a beating pattern in the Friedel oscillations, 
reveals many important details about the nature of interactions. Thus it has a potential to become 
an investigation tool in condensed matter physics. 

PACS numbers: 72.10.Fk, 72.25.Mk, 73.63.-b 



The role of electron-electron interactions is decisive in 
a large number of condensed matter phenomena. The 
correlations are not only responsible for such spectacular 
effects as the fractional quantum Hall effect and high- 
T c superconductivity, but they also turned out to influ- 
ence the physical properties of metallic ID and quasi-lD 
materials in a very profound way. In fact, contrary to 
their higher dimensional relatives, which are almost per- 
fectly described by the Fermi liquid theory, the ID metals 
even constitute a special universality class of Tomonaga- 
Luttinger liquids (TLL) [l], Q • Very naturally, it is ex- 
tremely important to have precise, reliable and easily 
accessible methods to investigate electronic interactions. 
The most useful of them involve quantities which are only 
present in interacting systems such as, e.g., the high en- 
ergy tails of field emission spectra 0, 0, @ • 

During the last two decades scanning tunneling mi- 
croscopy and atomic force microscopy as well as their 
descendants established themselves as some of the most 
efficient techniques for the investigation of low dimen- 
sional materials. It became possible to image the elec- 
tronic density of, e.g., carbon nanotubes, with atomic 
resolution In one of the recent experiments Q even 

an explicit imaging of Friedel oscillations (FOs) near the 
end of a single-wall carbon nanotube (SWNT) has been 
successfully achieved, making it possible to obtain the 
signatures of a TLL. 

The classical FO is nothing else but an electron den- 
sity modulation in vicinity of a hard wall potential Q. 
Due to the full reflection at a boundary at x — a co- 
herent superposition of incoming and outgoing electron 
wavefunctions results in a sm(2kp-x) / x behavior of the 
correction to the uniform electron density in the bulk (at 
x — > oo). The fundamental oscillation period then is de- 
fined by the Fermi wave vector kp and is equal to ir/kF- 
Similar density oscillations also emerge in the vicinity 



of any scatterer introduced into a metal. Their form and 
the actual particle density at the position of the scatterer 
of course depend on the scattering potential. One could 
conclude that in the case of a finite barrier transmission 
the modulation can change qualitatively so that a beating 
pattern due to two harmonic components would emerge 
as soon as the Fermi wave vectors on both sides of the 
scatterer become different, fcj, ^ fcjj. A straightforward 
calculation immediately shows that this is not the case as 
long as electron-electron interactions are not taken into 
account. Suppose we have a scatterer with reflection and 
transmission amplitudes r and t, respectively, placed at 
x = 0. Let dk and ak be the plane wave annihilation 
operators for the electrons injected into the system at 
x — > ±oo. Then the full field operator for electrons on 
the left side of the scatterer is 
dk 
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where for simplicity we assume the band structure en- 
coded in the density of states v(k) to be the same on both 
sides of the scatterer. The particle density is then given 
by the average {ij^ L {x)ip . Taking into account that 



cross correlations of the type 



4 dk> 



vanish, one real- 



izes that the term proportional to t*t is not ^-dependent. 
That is why the density oscillations on the (here) left side 
of the scatterer are not affected by the electrons from 
the right side up to an overall uniform x-independent 
shift. It is instructive to carry out the calculation for a 
specific system. We use the model for the conductance 
band of a SWNT: the electrons are assumed to be one- 
dimensional with a dispersion relation which is the same 
on both sides of the scatterer and is linear E(k) — vpk, 
where vf is the Fermi velocity (for simplicity we use units 
with h = e = m e = 1 throughout), see Fig. [T] Different 
values of hp on both sides can be produced by applying 
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FIG. 1: Sketch of the system under investigation. 

a finite bias voltage to the whole system complemented 
by appropriate gating of the constituent subsystems. For 
convenience we shall call the difference in the Fermi en- 
ergies 'voltage' V: Ef,l/R — v f^lir — vpkp ± V/2. 
Then the expectation values (a k ak') = ^-KSkk'n^k) and 
{df.dk') = 27r5fcfc'n^(fc) contain the Fermi distribution 
functions til/ r with Fermi energies Ef,l/r- The par- 
ticle density Pl/r is given by the sum of three terms: 
(i) the full constant average density p^ = kp/ir; (ii) a 
constant shift proportional to V on both sides given by 

p[ V /R = ±(1 - T)V/(2v f tt) , (2) 

the difference of which can be reworked into Lan- 
dauer's residual resistivity dipole strength p = kpV(l — 
T)/ (27m), where n is the density of the charge carriers 
(id 1 1 11 ] ; (iii) a space dependent oscillating part 

(osc), y/l-T sin(2x k L/R - yy) + sin(yy) 

Pl/r W - — ^ - • V) 

Here T = t*t is the transmission probability and ip r is the 
phase shift for a reflected electron. The full constant av- 



erage density p^ usually is compensated by the constant 
homogeneous background charge of the ions of the under- 
lying lattice structure and must therefore be subtracted 
in the interaction terms. Eqs. ([2|) and ((3]) show explic- 
itly that pl [pr] only depends on k^ [k R ] in the absence 
of interactions. Different ways to measure the induced 
chemical potential oscillations are in detail discussed in 

Now we want to address the effect of interactions on the 
FOs. We consider a spin-independent pointlike density- 
density interaction with strength U. In terms of the 
Keldysh Green's functions (GF) defined as g< L (x,y,t) = 
—i(t(} a (x,t)tjjL(y,0)} for the lesser GF and a similar def- 
inition of the retarded GF, the first-order correction to 
the particle density is given by 

8p L {x) =417 Im I dt ( dy (4) 

J J — OO 

x ^2 Pa{y) 9 r L^{x,ay,-t)g^ L {ay,x,t) , 

a=R,L 

with the prescription a = R,L = (±), which accounts 
for two different subsystems: the left one for y < and 
its counterpart on the right side, y > 0. This correction 
has a simple physical interpretation as it comes from the 
electron excursions from the measuring point x to some 
other coordinate y where it gets scattered on the local 
electron density p a (y) (which may be on either side of 
the barrier) and then propagates back to the original x. 
The scattering process is crucial as it allows for exchange 
of momenta. In this way the information about the Fermi 
wave vector on the other side of the impurity is carried 
back to the measuring point x. This effect is similar 
to secondary tunnelling processes (Auger like scattering) 
responsible for high-energy tails in field emission spectra 
of correlated hosts [H, 0, 11 • The GFs necessary for the 
evaluation of ((4]) are 



9 r L e l(x,y,t) = -iQ(-t) — e~ m - kt {cos[k(x - y)} - VT^T sm[k(x + y) - Vr }} , 

f°° dk f ~i 

9$ L {y,x,t) = \ nLe *Hy-*) + [ nL (i„ T ) + n R T]e^y-^ -2n L VT^ coslk(y + x)-<p r ]\e- 

Jo 2tti I J 



■ivp kt 



gT R (x,y,t) = -iQ(-t)— I dke- lVFkt sin[k(x + y) - <p r 



7T 







P OO 77 

g^ L (y,x,t) = - J — [n L VTe'^y +x ^ - n R VTe^y +x ^ - y/T(l - T) (n L -n R ) e ik ^} e 



—ivpkt 



i 

where 9 denotes the unit step function. Note that these the exact scattering states. As we are looking for the sec- 
GFs are exact in scattering strength since we are using ondary FO (contrary to the primary one coming from the 
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native fcj,) it is more convenient to perform the Fourier 
transform of the correction to the density. Furthermore, 
we would like to concentrate on the consequences of hav- 
ing two different Fermi levels in the system, thus we only 
keep terms proportional to the difference of the Fermi 
distribution functions An = (nj, — tir). 

The detailed evaluation of ([4]) reveals two contribu- 
tions to the first order density correction 5 p. The first 
one comes from the mean shift p^) R of the zeroth order 
density due to a finite V and it is found to be given by 



s?ZHp) = T(1 - T)VU 



4(irv F ) 2 
1 + e ~ 2i ^ 



(5) 



■jrQ{& + 2k L )Q(-&-2k R ) 
2,k n 



-z£]T alog 1 + ^ [f-e(£p + fc Q )(l + e 2i ^)]) 

a=R,L ^ P ' 

The spectrum of the density profile shows up two distinct 
peaks at p = ±2fc£/^j while in the noninteracting case 
there would be only one native periodicity with wave- 
length 7t//cl . The peaks are in fact logarithmic singular- 
ities due to the one-dimensional geometry of the system. 
By construction the signal strength is proportional to 
the applied voltage. It vanishes for both perfect trans- 
mission, when translational invariance of the system is 
restored and no buildup of FOs is possible, and perfect 
reflection, when the electrons on the one side of the bar- 
rier cannot acquire information about the Fermi edge on 
the other side. The prefactor T(l — T) can indeed also 
be explained more intuitively. Suppose we have a wave 
packet approaching the barrier coming from the left reser- 
voir. At the boundary a part of it is reflected with an 
amplitude ~ y/l — T. The transmitted part which is 
~ s/T propagates to the point y where it interacts with 
the wave packet coming from the right reservoir, which 
on its way to y is reflected from the boundary picking up 
a prefactor ~ y/1 — T. The interaction process yields a 
prefactor U and the propagation back through the bar- 
rier gives another VT, see Fig. [2j Thus the full prefactor 
is (VTVl - T) 2 = T(l - T). For a representative set of 
parameters the density spectra as given by ([5]) are plotted 
in the left panels of Fig. [3] Being transformed back into 
the coordinate space, the two peaks result in a beating 
pattern of the density oscillation as a function of distance 
from the impurity. 

Besides these oscillations, originating from the finite 
(average) coordinate-independent difference in chemical 
potentials, there is, of course, a contribution from the 
oscillating part of the density on the other side of the 
barrier. The physical motivation and qualitative depen- 
dence on the transmission is the same as discussed above, 
but corresponding expressions are rather long, so we re- 
strict ourselves to the plots shown in the right panels of 
Fig- 131 As before, distinct features arise around both dou- 




FIG. 2: Graphical representation of the processes contribut- 
ing to the secondary Friedel oscillations, see text. 
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FIG. 3: Real and imaginary parts of the density oscillations 
spectrum (arbitrary units) for V = O.lfc^ 1 . Left panels: con- 
tributions from Right panels: contribution from the 
oscillating part of the density pL° ■ 



bled Fermi momenta p = ±2k L / R . We expect that this 
interaction signature should be observable in dedicated 
experiments, thereby giving important insights into the 
nature of electronic correlations. 

So far we treated the system exactly in the tunnelling 
(scattering) amplitude and performed a weak interac- 
tion expansion. In strictly ID environments, at least 
in equilibrium, the electron system belongs to the uni- 
versality class of TLLs. In this case the interactions 
have to be treated nonperturbatively. For the FOs at 
an open boundary (at a hard wall) the density profile 
can be calculated analytically and numerically for arbi- 
trary interaction strength [l3|> [3, EH- In our case the 
situation is much more difficult as we must allow for at 
least weak electron tunnelling at a barrier between two 
half-infinite TLLs. It turns out, that such a situation in- 
deed describes the low-energy fixed point of a TLL with 
an impurity li|, 17 1 . Therefore we adopt this model and 
perform a calculation which is 'dual' to the one presented 
above, taking into account the interactions exactly while 
performing a perturbation expansion in the tunnelling 
amplitude. 

The beating pattern in the density oscillations already 
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shows up in the inhomogeneous case (which is more con- 
venient to handle), when a noninteracting subsystem is 
coupled to a genuine TLL and the FOs are measured in 
the interacting part The correction to the density profile 
is then found to be given by 



6p(x) 



M2k F + V/v F )x g ^ (^)) . ( 6 ) 



The evaluation of the expectation value on the rhs is 
accomplished with the help of the bosonization represen- 
tation of the TLL Hamiltonian with the interaction of 
the kind we already have used in the first part of this 
Letter. The correlation strength translates into the di- 
mensionless interaction parameter g = 1/ Wl + U/(ttvf), 
with < g < 1 for repulsive interaction. The correction 
to the density is 



8{1>l{x)M-*)) = ~ 



-1/9 



ttvf (ttoo) 2 1 2a3 |s 
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where ao is the lattice spacing of the underlying lattice 
model. F(x) is an oscillatory function, whose Fourier 
transform can be written in the form 




+ T-/2 , fe , 
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As 1/g > 1 for repulsive interactions it is the in- 
dependent parts in the brackets which dominate the be- 
haviour of F(k). A stationary phase argument with re- 
spect to the r_ integration yields the condition k/2 + 
V/vf = for the maximum of F(k). This is confirmed 
by the numerical evaluation of |7J) shown in Fig. SI where 
F(k) is plotted for different values of the interaction pa- 
rameter g. Therefore we conclude that the correction 
(|rj|) shows up oscillations with period it/k^ in addition 
to the native one with ir/kr, just like in the perturbative 
interaction calculation. Interestingly, the stronger the in- 
teractions the more pronounced are these newly arising 
oscillations. Similar effects in a slightly different setup 
have been reported in Ref . [HI . 

To conclude, we have discussed the density profile of 
Friedel oscillations in a ID metal with different Fermi 
wave vectors on either side of an imbedded impurity. 
Contrary to the naive expectation in the absence of inter- 
actions, the density oscillations possess only one period- 
icity originating from the local value of kp even when the 
impurity has a finite transmission coefficient. The situa- 
tion changes dramatically in interacting systems, where 
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FIG. 4: The function F[k) as given by Eq. © plotted for 
different values of the interaction parameter g. 



the spectral function acquires an additional peak, which 
is a result of the interplay of interactions and the Lan- 
dauer dipole. We have shown this in two different situa- 
tions by first taking into account interactions perturba- 
tively and using the exact scattering states. In a second 
step we have considered the dual situation of a weak link 
between two TLLs where the electronic correlations have 
been taken into account with the help of bosonization 
technique. We are convinced that the very fact, that the 
beating pattern owes its existence to the interactions has 
the potential to make measurements of Friedel oscilla- 
tions (which by now have become unproblematic) an im- 
portant tool to study electronic correlations in ID metals 
and carbon nanotubes in future experiments. 
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